Gauge Fields, Geometric Phases, and Quantum Adiabatic Pumps 
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Quantum adiabatic pumping of charge and spin between two reservoirs (leads) has recently been 
demonstrated in nanoscale electronic devices. Pumping occurs when system parameters are varied 
in a cyclic manner and sufficiently slowly that the quantum system always remains in its ground 
state. We show that quantum pumping has a natural geometric representation in terms of gauge 
fields (both Abelian and non-Abelian) defined on the space of system parameters. We make explicit 
the similarities and differences with Berry's geometric phase. Tunneling from a scanning tunneling 
microscope tip through a magnetic atom could be used to demonstrate the non-Abelian character 
of the gauge field. 
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Normally transport of electrical charge is dissipative 
(i.e., it produces heat). However, quantum adiabatic 
pumping [lj provides a means in nanoscale electronic de- 
vices to use novel quantum effects to transport single 
electrons with minimal dissipation [2(. Furthermore, it 
is also possible to pump electron spin without pumping 
charge 0, 3 . Both charge Q and spin pumping have 
been recently achieved experimentally, by cyclic variation 
of the gate voltages that control the shape of an open 
quantum dot. This motivated extensive theoretical re- 
search in this topic, especially on quantum charge pump- 
ing |?l IM 1^. llOl Quantum spin pumping opens the 
way for applications in spintronics. It is sometimes sug- 
gested, but not explicitly shown, that quantum pumping 
is related to Berry's phase. As first emphasized by Berry 
|12| , discrete quantum systems have the counter-intuitive 
property that when some of the parameters controlling 
the system are slowly varied and brought back to their 
initial value the quantum state of the system is different 
to the initial state. That is, a quantum state may acquire 
a geometric phase exp(i7c) in addition to the normal dy- 
namic phase exp(— (i/h) J E(t)dt). This remarkable dis- 
covery may be recast in the language of holonomy theory 
[l3j| . Subsequent work showed that non-Abelian gauge 
potentials can arise as a result of degeneracies of energy 
levels of the system [ill 111 Il6| . However, it is unlikely 
that adiabatic pumping, characteristic of quantum open 
systems considered in this Letter, results from Berry's 
phase for closed discrete systems. 

We present a systematic treatment of quantum adia- 
batic pumping in open systems in terms of parallel trans- 
port and gauge fields (both Abelian and non-Abelian) 
defined on the system parameter space, which reveals a 
unifying concept of geometric phase underlying scatter- 
ing states. We make explicit the similarities and differ- 
ences with Berry's phase associated with cyclic variations 
of closed quantum systems (both degenerate and non- 
degenerate) (see Table 1). In the scattering approach 
developed by Brouwer |17j . based on an earlier work of 
Biittiker, Thomas, and Pretre [l^, a compact formula 
was presented for the pumped charge (current) in terms 



of the parametric derivatives of the time-dependent scat- 
tering matrix subjected to the modulating potential. We 
show that the pumped charge, given by Brouwer's for- 
mula 17], is essentially the geometric phase associated 
with the U(l) subgroup of the gauge group U(M) (M 
is the number of channels in a certain lead), whereas 
the non-Abelian sector SU(M) describes the adiabatic 
pumping associated with the internal degrees of freedom 
such as spin. Expressions are given for the gauge poten- 
tials associated with tunneling from an STM (scanning 
tunneling microscope) through a magnetic atom. We 
suggest an experiment which can be used to illustrate 
the non-Abelian character of the gauge field. 

The quantum system. Consider a mesoscopic system 
with N leads, and for the n-th lead there are M n chan- 
nels. Our aim is to study quantum pumping by period- 
ically varying a set of the independent external param- 
eters X = (X , • • • , X v , • • • , X p ) slowly as a function of 
time t. In the scattering approach, the S matrix is an 
JV x TV matrix with M the total number of channels, 
N = J2n=i We define vectors n a = (S a i, •• ■ , S a _\f) 
in terms of the rows of the scattering matrix S[X(t)] 
associated with the n-th lead. The unitarity of the scat- 
tering matrix implies that these vectors are orthonormal 
n* • = 5 a p, a,f3= 1, • • • , M n . That is, this provides 
us with a smooth set of (local frame) bases n a (t). 

The gauge potential. Assume that the parallel trans- 
port law 











(1) 



holds for (fibre) vectors ^f a , where d^ a is the variation in 
resulting from a variation dX in the external param- 
eters. If \1/q(0) = n Q (0), i.e., the initial vector describ- 
ing the scattering process in which the incident particle 
comes from the cv-th channel in the n-th lead, then the 
degeneracy of the channels implies that the transported 
vector ty a (t) must be a linear combination of all n a (t), 
^a(t) = TlipU a i3{t)r\.p(t). Expressed another way, the 
transported vector describes a combined scattering pro- 
cess in which particles come from all channels in the n-th 
lead. Obviously, U(t) is unitary, i.e., U(t) £ U(M n ). 
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Physically, this means certain information about where 
the incident particles come from is lost during parallel 
transport, and is encoded in the unitary matrix U(t). 
Inserting into the parallel transport law in Eq. ifTJl. we 
have 

(U- 1 dU) aP = -n} ■ dn a . (2) 

Since n Q varies as the parameters X u vary with time, we 
can thus define the gauge potential A a p v = • <9„n Q , 
where d v = djdX v so that 

(U- x dU) af3 = -J2A a0 „dX» . 

V 

This can be integrated in terms of exponential integrals. 
For the period r of an adiabatic cycle, we have 

(C/(t)) q /3 = (P exp(- I J2 A v dX v )) a0 , (3) 

where P denotes path ordering. Defining A = 
^2 V A v dX v , one can see it is Lie algebra u{M n ) valued 
and thus anti-Hermitian. A plays the role of a gauge 
potential, as in the case of Berry's phase 0] for closed 
(discrete) quantum systems. 

Gauge transformation. The gauge group U{M n ) orig- 
inates from the unitary freedom in choosing local bases 
n a (a = 1, ••• ,M„), 



This amounts to different choices of the scattering ma- 
trix: S'(t) = fl(t)S(t) with fi(t) a diagonal block ma- 
trix, the n-th block of which is an M n x M n unitary 
matrix uo(t). Physically speaking, left multiplication of 
the scattering matrix S(t) by £l(t) just redistributes the 
scattering particles among different incoming channels 
associated with a certain lead, which does not affect cor- 
relations at the scatterer and so the physics remains the 
same. The gauge potential A(t) transforms as 

A'(t) = dujuo- 1 + ujAuj- 1 . (4) 

The gauge field strength defined by F = dA— AAA trans- 
forms covariantly F' = loFlo -1 . Therefore, a U(M n ) = 
U(l) x SU(M n ) gauge field is defined on a p-dimensional 
parameter space which drives the quantum pumping. 
The trace of U(t) given by Eq. J3J is gauge- invariant. 

Justification of parallel transport. Now we need to jus- 
tify the assumption of the parallel transport law in Eq. 

Physically, by "adiabatic" we mean that the dwell 
time t p during which particles scatter off the scatterer 
is much shorter than the time period r = 2tt /uj a during 
which the system completes the adiabatic cycle. Here Lo a 
is a slow frequency characterizing the adiabaticity and 
t p is related to (but not determined alone by) Wigner 
time delay matrix t w (s,E) = —iS^ (s, E)dS(s, E)/dE 



with E being the energy of scattering particles and s be- 
ing the so called epoch defined as s = u> a t [lflj . Then 
the response to the variation of the particle distribu- 
tion in a certain channel is only limited to channels as- 
sociated with the same lead. That means we ignore 
any responses which involve channels in different leads. 
Such a response can be treated as dissipation, a correc- 
tion to the adiabatic limit. Then, the parallel trans- 
port law in Eq. follows from the parallel transport 
law for the wave function. The latter is a solution of 
the time-dependent Schrodinger equation ihdt\<$>(t)} — 
H{t)\<fr(t)) . As is well known, the Schrodinger equa- 
tion induces a parallel transport law lm((p(t)\dt(j>(t)) = 
20], with = exp(i / h(t)dt)\$(t)), where h(t) = 

We can write the wave function \(j)(t)} 
as a linear combination of all scattering states associ- 
ated with a certain lead in the adiabatic case. Formally, 
l<K*)) — Eq c a\^a{t))i with \ip a (t)) denoting the scat- 
tering states in which the scattered particles come from 
channels associated with the n-th lead, and c a being ar- 
bitrary constants. Then we have (tp a (t)\dtiljp(t)) = 0. 
The adiabatic assumption implies that \ipa(t)} may be 
expanded in terms of instantaneous asymptotic scatter- 
ing states, |^ a (t)> = ^ U af3 (t)\^ (t)), with |V>f(t)) = 

\ot)in + Y$f=i Sap(t)\P)out, a = 1, • • • , M n . Here, \a) in 
and 1/3} o«t denote, respectively, the incoming and outgo- 
ing scattering states, which are normalized such that they 
carry a unit flux. Substituting into the parallel transport 
law for \ip a ), one gets Eq. J5J which is equivalent to the 
parallel transport law for row vectors of the scattering 
matrix. 

Quantum adiabatic pumping. In order to establish 
the connection between the geometric phase above and 
the quantum pumping charge, we need to consider 
the time-reversed scattering states (t)) — \a) in + 
Y^=i Sp a (t)\/3) ou t with * denoting the counterparts un- 
der time reversal operation |2l| . which constitute a so- 
lution of the Schrodinger equation for the time-reversed 
Hamiltonian H at any given (frozen) time at the epoch 
scale 19]. This gives rise to another gauge potential 
A a p v = • d u h a with n Q = (S la , ■ ■ ■ ,Sx a ), i-e., the 
column vectors of the scattering matrix S(t). In this 
case, the gauge group arises from redistribution of scat- 
tering particles among different outgoing channels. If 
= S/3 w a /9n^(t), then the gauge transformation 

takes A'(t) = duitli^ 1 + ujAuj^ 1 . The gauge fields A 
and A are connected via time reversal operation. If we 
identify the emissivity into the a-th channel in the n-th 
lead as Im[A aa /27r]_L8|, then we immediately reproduce 
Brouwer's formula [17| describing charge pumping, which 
turns out to be associated with the Abelian subgroup 

Q=^-Im^Tri, (5) 

with Q being the charge transferred into the n-th lead 
during one cycle. That is, the charge transferred dur- 
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ing adiabatic pumping is essentially the geometric phase 
associated with the charge sector U(l). This also ex- 
plains why Planck's constant H does not occur in the 
adiabatic quantum pumped charge (current), a pecu- 
liar feature different from the Landauer-Biittiker con- 
ductance. However, as is well known, the geometric 
phase is determined only up to a multiple of 2ir. This 
concerns global geometric properties, i.e., the winding 
number of the overall phase of the gauge transformation 
in Eq. Q), N = l/(2ni) § Tr^wuT 1 ). The require- 
ment that all physical observables be invariant under 
the gauge transformation leads us to the conclusion that 
0(Q) = 0{Q — eN), with O denoting any observable. 
This result has been noticed by Makhlin and Mirlin , 
without proper justification, for the counting statistics in 
quantum charge pumps (see also, Ref. 0). 

To see the effects caused by non-Abelian gauge po- 
tentials, we need to consider gauge invariant quanti- 
ties. Besides TrC/(r), we see that both the determinant 
and eigenvalues of U(t) are gauge invariant. Actually, 
there are M n independent gauge invariant quantities such 
as the eigenvalues exp[i7 Q ]. On the other hand, there 
are M n independent simultaneous observables such as 
the pumping currents I a flowing into the a-th channel, 
which must be gauge invariant. Therefore, one may ex- 
pect that the pumping currents I a are some functions 
of 7i,-- - , 7m„- Because our argument only relies on 
gauge invariance and does not depend on any details of 
the system, such functions must be model-independent. 
Guided by the results for the so-called "Abelianized" 
non-Abelian gauge potentials, i.e., the gauge potentials 
which turn out to be diagonal in a certain gauge, we have 

I a = _ 2^rT«- ( 6 ) 

This connects the physical observables with the eigen- 
values of the geometric matrix phase. Especially, the 
charge pumping current I c corresponding to the Abelian 
sector is I c = Y, a I a = -l/(27tr) Y, a 7a- One 

may verify that this is consistent with Eq. © since 
Q(t) = el c and ddetU(t) = -Tri(i) det U(t). Alter- 
natively, Q — e Im In det U (r) . Similarly, we may define 
generalized "spin" pumping currents associated with the 
Cartan subalgebra of the non-Abelian sector SU(M n ). 
The simplest non-Abelian case U(2) is relevant to the 
charge and spin pumping. 

Tunneling through a single magnetic spin. Consider 
the Hamiltonian which describes a system consisting of 
two leads coupled to a single site, the spin of which has 
an exchange interaction J with a magnetic spin |22j , 

+ £ (V4,4'+H.c). (7) 

Here cl and Cka are, respectively, the creation and de- 
struction operators of an electron with momentum k and 



spin a in either the left (L) or the right (R) lead, and d\ 
and d a are the counterparts of the single electron with 
spin a at the spin site. The quantity e k a are the single 
particle energies of conduction electrons in the two leads, 
which we will assume tka = Vp{\k\ — Hf) with the conven- 
tion that Vf = li and the momentum is measured from 
the Fermi surface for electrons in leads. The electrons 
on the spin site are connected to those in the two leads 
with the tunneling matrix elements Vu a ,a' ■ For simplic- 
ity, we assume symmetric tunneling barriers between the 
local spin and the leads, and only keep the spin-conserved 

coupling; viz. Vl++ = Vl = Vr++ — Vr— = V and 

V L +- = V L -+ = V n +- = Vr-+ = 0. The entries of the 

coupling matrix Q take the form = — f2 = cos# 

and f2_| = Q^_ + = sin#exp(— i(j)). The model is exactly 

soluble as far as the scattering matrix is concerned. 

Once the scattering matrix is determined, our general 
formalism leads us to the non-Abelian gauge potential, 

A = A e d8 + A^dcj), (8) 

where Ag = Ala 1 + Aja 2 + A^a 3 with A\ = i(sin(£i - 
82) cos cos (/)+ (1 — cos(#i — S2)) sin 0)/4, A 2 , = i(sm(5i — 
62) cos8sin(j) — (1 — cos($i — <5 2 )) cos(/>)/4, and Ag = 
-i(sin(<5i - 8 2 ) sin 6»)/4, and i = A\o x + A\a 2 + i^cr 3 
with A]p = — i(s'm(Si — 62 ) sin 8 sin <p — (1 — cos(<5i — 
£2)) sin^cos0cos^)/4, A\ — i(sin(<5i — <5 2 ) sin0cos(/) + 
(1 — cos(#i — 82)) sin 8 cos 8 sine/)) /4, and A^ = — i(l — 
cos(<5i — 82)) sin 2 9/ A. Here Si (i = 1,2) are the phase 
shifts defined by 81 = — 2tan~ 1 (r/(fc — J)) and 82 = 
-2tan" 1 (r/(fc + J)) with the tunneling rate T = V 2 . 

The gauge field strength F then takes the form 
F = - cos(<5i - 5 2 ))n ■ adQ/A. Here dVt = 

sin 8dd A d<j) is the invariant area element and n — 
(sin 8 cos 4>, sin 8 sin cf>, cos 8) is the direction of the mag- 
netic spin. Obviously, this is just a simple rotation of 
the standard form F = —i(a 2 — l)a 3 dil/2. Up to a 
gauge transformation, this is the same non-Abelian gauge 
potential, found by Moody et al. 15] for a diatomic 
molecule. This is consistent with a theorem, proved in 
|23| , stating that the rotationally invariant connection on 
the sphere is essentially unique. To establish the relation 
between a and cos(<5i — ^2), we need to calculate the gauge 
invariant quantity Tri^ A with *F being the dual of 
F. Then we have a 2 = (3-cos(<Si -5 2 ))/2. When J = 0, 
i.e., in the absence of the direct exchange interaction be- 
tween electrons and the local spin, the gauge field is a 
pure gauge because F = 0. Since the Pauli matrices are 
traceless, we have Tr^4 = 0, meaning that charge pump- 
ing is absent in the model under consideration. That 
implies 7+ = —7-. Therefore, the spin pumping current 
defined by I s = 1+ — I- becomes I s = —j+/(ttt). 

One can compute a phase factor Usr. which is ob- 
tained from the time-reversed counterpart of Eq. J3J for a 
"spherical rectangle (SR)". From \/2TcUsr — cosittI s , 
the spin pumping current / s may be extracted and is 
shown in Fig. ^ as a function of system parameters for 



4 



two different paths, C\ and C%. The spin pumping cur- 
rent takes its maximum value around the resonant scat- 
tering lines k = ± J. 

Possible experiments. A scanning tunneling micro- 
scope (STM) has been used to detect a quantum mi- 
rage around a single magnetic cobalt atom placed on a 
non-magnetic metallic copper surface |24l . Electron spin 
resonance (ESR)-STM experiments [25[ have advanced 
to the point that they have spatial resolution at the level 
of a few spins j2||. The STM setup as shown in Fig. 
should make it possible to observe gauge invariant spin 



pumping via a single magnetic atom on the surface of the 
substrate. To measure the spin pumping current / s , one 
could replace one of the leads by a ferromagnetic one. 
The spin pumping current can then be measured via the 
charge pumping currents. 
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Berry's Phase 


Scattering (Pumping) Geometric Phase 


Closed systems 


Open systems 


Wave functions \ip a ) 


Row (column) vectors rio, (n a ) of the S matrix 


Energy levels E n 


Leads n 


M degeneracies 


M channels 


Discrete spectrum (bound states) 


Continuous spectrum (scattering states) 


Eigenstates in the n-th level 


Channels in the n-th lead 


Parallel transport due to adiabatic theorem 


Parallel transport due to adiabatic charge and spin pumps 


rjauffp nntpntial A~a,. — ( ?/;^ r),, \ih~ > 


Gauge potentials A a ^^ — • Di/Yia and A a p u — fi^ ■ O^Yia 


U(M) bundle 


U(M) bundle 


Gauge group U(M) arising from 


Gauge group U (M) arising from redistribution of the 


different choices of bases 


scattering particles among different channels 


External parameters X — (X 1 , • • ■ , X p ) 


External parameters X — (X 1 , • ■ • , X v ) 



TABLE I: Comparison of Berry's phase and the quantum scattering (pumping) geometric phase. 
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FIG. 1: The dependence of the spin pumping current I s (times 
r, the period of cyclic variation) on system parameters for 
two leads connected to a single magnetic spin whose direc- 
tion is slowly varied around the path shown on the right. 
A. Left: Schematic of the magnetic spin coupled to left (L) 
and right (R) leads. The magnetic spin S precesses around 
the direction of the magnetic field B. Right: An equiva- 
lent scanning tunneling microscope experimental setup. The 
pumping cycles on the parameter (9, 0)-sphere are, respec- 
tively, taken to be Ci(#i, 01, 02, t/>2) = (t/8, 0, tt/2, it) for B 
and C 2 (9i,4>i,e 2 ,(t>2) = (t,/8, 0, 7t,/8, tt) for C. 



